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Abstract. We show that a hydrodynamical model with continuous particle emission 
instead of sudden freeze out may explain both the observed strange particle and pion 
abundances and transverse mass spectra for light projectile at SPS energy. We found 
that the observed enhancement of pion production corresponds, within the context of 
continuous emission, to the maximal entropy production. 
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1. Introduction 

The main purpose of tlie ongoing and future heavy ion programs at the high energy 
laboratories (CERN, BNL) is to investigate the formation and properties of hot dense 
matter, in particular, the phase transition from hadronic matter to quark gluon plasma 
(hereafter QGP) predicted by Quantum Chromodynamics (QCD). Various possible 
signatures of the appearance of the QGP have been suggested: entropy increase (due 
to the release of new degrees of freedom, namely color), strangeness increase (due to 
enhanced strange quark production and faster equilibration), J/ip suppression (due to 
color screening or collision with hard gluons) or enhancement at higher energy (due to 
recombination of dissociated cc pairs), production of leptons and photons (emitted from 
a thermalized QGP and unaffected by strong interactions), etc. These signals have been 
studied extensively in experiments (see for example PJ). More recently, the observed 
high transverse momentum depletion in Au+Au central collisions at RHIC energies 
(jet quenching) and also its absence in d+Au collisions, together with the system-size 
dependence of mono-jet formation, have been considered as a convincing evidence of 
the formation of QGP at RHIC 2J. On the other hand, some of the above mentioned 
signals are purely of thermodynamical nature and sometimes their significance is not well 
defined for finite systems. In such cases, a more careful analysis would be necessary to 
clarify the effect of finiteness and dynamical evolution of the system on these signals [3]. 
Some authors suggest that the incident energy dependence of several quantities should 
be studied and they claim that the set of incident energy dependences of particle 
multiplicity, average transverse momentum and kaon-pion ratio as a whole indicates 
the appearance of the mixed phase in central collisions at SPS energies j3]. In fact, 
if we consider one specific signal of thermodynamical nature, there are many factors 
which may give a similar response of the claimed signal. For example, it is well known 
that the hadronic final state interaction also works to suppress J/ip and the system 
size dependence should be carefully studied to extract the significance of the observed 
data. Therefore, it is always important to study and grasp the hadronic effects on other 
proposed signals. 

A major problem to trace back any signature unambiguously to a quark gluon 
phase is that it is still unknown which theoretical description describes best high energy 
nuclear collisions. On one extreme, one might use a microscopic model. In principle, 
such approach would provide a faithful realization of the true physical processes if all the 
relevant physical degrees of freedom and corresponding interactions were incorporated 
appropriately. However, in practice, this is not possible and some hypothesis and 
simplifications are necessarily introduced. For example, we may restrict ourselves to 
pure hadronic degrees of freedom, but is is known that such models fail[6j to reproduce 
simultaneously strange and non-strange particle data in nucleon-nucleon collisions and 
central nucleus-nucleus collisions at SPS energies|. Partonic microscopic models are 

I Modifications have been attempted to solve the strangeness problems in hadronic microscopical 
models, see however [3] 
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expected to work at energies higher than SPS§. One promising approach to describe the 
final state interactions in heavy ion colhsions is the UrQMD model jS] . However, due to 
the complexity of the calculation and many uncertainties of input data such as hadronic 
cross sections, it is not easy to get a simple physical insight for the behaviour of the 
observed quantities. 

On the other extreme, one might use a thermal or hydrodynamical model. In 
such models, it is assumed that a fireball (region filled with dense hadronic matter or 
QGP in local thermal and chemical equilibrium) is formed in a high energy heavy ion 
collision and evolves. Hydrodynamical models have been used successfully to describe 
various kinds of data from AGS to RHIC energies. In particular, at SPS, they are 
able to account for strangeness data but, in some simplest versions, fail to predict large 
enough pion abundances (cf. next section). However such an approach has a great 
advantage compared to the microscopic models in the sense that very few input data 
are necessary because of the assumption of the local thermal equilibrium. Furthermore, 
we know that this macroscopic description works quite well suggesting that the local 
thermal equilibrium in heavy ion collisions is a reasonable approximation. Thus one 
might ask whether it is possible to introduce a correction to the local thermodynamical 
equilibrium and its effects on the observed particle abundances. The aim of this paper 
is to study this problem. We will argue that, within a very simple and natural picture of 
particle emission (continuous emission), the observed particle abundances and spectra, 
including strange particles are consistently reproduced. In addition, we show that 
the non-equilibrium component leads to an increase of entropy of the system, hence 
a substantial pion increase compared to the usual sudden freezeout mechanism. 

2. Hydrodynamical description with (standard) freeze out emission 

In the standard hydrodynamical models, one assumes that particle emission, called in 
this case freeze out, occurs on a sharp three-dimensional surface (defined for example by 
T{x,y, z,t) = constant). Before crossing it, particles have a hydrodynamical behavior, 
and after, they free-stream toward the detectors, keeping memory of the conditions 
(flow, temperature) of where and when they crossed the three dimensional surface. The 
Cooper-Frye formula jH] gives the invariant momentum distribution in this case 



daf^ is the surface element 4-vector of the freeze out surface a and / the distribution 
function of the type of particles considered. This is the formula implicitly used in all 
standard thermal and hydrodynamical model calculations. It can be integrated to get 
the particle abundance, which depends only on the freeze out parameters ||, e.g. Tj out- 

§ However see and references therein. 

11 In the past few years, models assuming two different freeze outs, respectively chemical (abundances 
fixed) and thermal (shape of spectra fixed) have been frequently used. (It is then necessary to modify 
JH) cf. [ini.) At RHIC, this view is being challenged |TT]. 




(1) 
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Freeze out parameters can be extracted by analyzing experimental particle 
abundances. This has been done by many groups (for a review see e.g. [T^). The 
models have some variations among them in particular some fit 47r quantities while 
others consider quantities in fixed rapidity window, each approach having its own 
quahties and draw-backs. For a number of these approaches, it was noted that while 
they can reproduce strange particle abundances, they underpredict the pion abundance. 
This was first noted by ^] in a study of NA35 data and emphasized by O ^] 
in an analysis of the WA85 strange particle ratios and EMU05 specific net charge 
Dg = (N^ — N^)/(N~^ + N~) (with A^"*" and N~, the positive and negative charge 
multiplicity respectively). A similar problem arises with the Pb+Pb data from NA49 

[mini- 
Various possible improvements have been suggested so that these models could 
yield both the correct strange particle and pion multiplicities: sequential freeze out [IHj, 
hadronic equation of state with excluded volume corrections 1201 HI] non-zero pion 
chemical potential [131 1201 HOj, equilibrated plasma undergoing sudden hadronization 
and immediate decoupling [T3J [12112115 etc. In this paper we follow a different strategy. 
We feel that the assumption of sudden freeze out on a 3-dimensional surface is a drastic 
one; in addition it is not sustained by simulations using microscopic models [221- So we 
study a different particle emission mechanism, continuous emission. 

Before we turn to this, let us note that, to compare particle abundances in the 
continuous emission and freeze out scenarios, we will use a simplified framework to 
describe the fluid expansion, namely we suppose longitudinal expansion only and 
longitudinal boost invariance [221 • This approximation allows to carry out some 
calculations analytically and turns the physics involved more transparent. It is implicit 
however that this description applies at best to the midrapidity region and light 
projectiles (for S+S data, transverse expansion must be small, see e.g. [211 )■ should 
therefore consider midrapidity data such as S+S data from NA35 or WA94, or S+W 
data from WA85. In fact since we want to consider strange particles and non-strange as 
well, we will concentrate on NA35 data. In this simplified framework, in the case of a 
fluid with freeze out at a constant temperature and chemical potential, the Cooper-Frye 
formula (0) can be re-written ignoring transverse expansion as [2^1 

r- = ^rf,outiTf,out,To,To)m^}^iT) exp(— )i^i(- )(2) 

ayp±ap±ly=o n=l J-f.out J-f.out 

(The plus sign corresponds to bosons and minus, to fermions.) It depends on the 
conditions at freeze out: Tf ^ut and /i/.out = /^t f.outB + fis f.outS, with B and S the baryon 
number and strangeness of the hadron species considered, and fJ^s f.out{f'bf.out,Tf_out) 
obtained by imposing strangeness neutrality. So the experimental spectra of particles 
teach us in that case what the conditions were at freeze out. 
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3. Hydrodynamical description with continuous emission 

The notion that particle emission does not necessarily occur on a three dimensional 
surface but may be continuous was incorporated in a hydrodynamical description in 
[26] . In this model, the fluid is assumed to have two components, a free part plus an 
interacting part and its distribution function reads 

fix,p) = ffreeix,p) + fint{.X,p). (3) 

ffree couuts all the particles that last scattered earlier at some point and are at time 
in X. fint describes all the particles that are still interacting (i.e. that will suffer 
collisions at time > x"). The invariant momentum distribution is then 

Ed^N/dp^ = J d^xD^[p^ffree{x,p)]. (4) 

D ^j\p^ f free{x 1 p)] is a covariaut divergence in general coordinates and d^x is the invariant 
volume element. A priori formula (j3)) is sensitive to the whole fluid history and not just 
to freeze out conditions as in formula (^. 

In the simplifled framework that we will use to compare particle abundances in 
the continuous emission and freeze out scenarios, we can approximate the equation of 
continuous emission (jlj) as 

dN 2g r m± cosh rjTp pdp + p j_ cos (ppprdr 

dyp±dp±^y^Q (27r)2 j-p=o.5 exp((m_L cosh 77 — /i)/T) ± 1 

where V is the probability to escape without collision calculated with a Glauber formula, 
Tp and pf are respectively solution of V{tp, p, 0, rj; v±) = 0.5 and V{t, pp, 0, 77; v±) = 0.5 
where v± is the particle transverse velocity. In (0), various T and p = p^B + psS appear 
{ps is obtained from strangeness neutrality), reflecting the whole fluid history, not just 
Tf out and pb f,out- This history is known by solving the hydrodynamical equations of 
a hadronic gas with continuous emission; it depends only on the initial conditions Tq 
and pu) at the initial time tq (we use the standard value Tq = 1 fm) . Therefore (0) 
only depends on the initial conditions. We expect that heavy particles, due to thermal 
suppression, will bring information on early times when they are more numerous. Fast 
particles can escape more easily from dense matter, so they will probe early times, 
whereas slow particles will probe late times, when they are flnally in diluted matter and 
make their last collision. Light particles such as the pion should probe the whole fluid 
history. 

4. Results for continuous emission 

In ^22j, we showed that it is possible to flnd initial conditions for the hydrodynamical 
evolution of the fluid that leads to WA85 strange particle ratios. In this paper, we want 
to broader this study. First we want to show that not just ratios, but abundances can 
be reproduced for a certain set of initial conditions. This is not trivial: for freeze out, 
fitting abundances or ratios is not very different because of volume cancellation, but 
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not for the continuous emission which is a process sensitive to the whole fluid history, 
particle mass and particle velocity. Second, we want to check that for the same initial 
conditions that reproduce abundances, spectra can be reproduced. 

Abundances can be obtained by integrating (5). Figure 1 shows the allowed 
region of initial conditions that lead to the experimental NA35 y = values |2H! 

A = 1.26 ± 0.22p^>o.5GeV', A = 0.44 ± 0.16p^>o.5GeV', Kl = 1.30 ± 0.22p^>o.62GeV, 

/i^ = 26 ± 1 and p — p = 3.2 ± 1.0. We do not use the and K~ abundances 
because they were measured outside the mid-rapidity region. Considering only strange 
particles, the allowed window is Tq ~ 183 — 188 MeV, /if,o ~ 70 — 125 MeV, with an 
ideal hadron gas equation of state and the strangeness saturation factor (a multiplicative 
factor for (5)) 7^ = 1.3. Including the p — p abundance decreases the window for /i^o 
to ~ 70 — 105 MeV. Finally including the abundance of negative particles leads to the 
small window Tq ~ 185 MeV and fibo ~ 100 MeV. 

Using a more soflsticated equation of state, the value of Tq might be decreased 
by some 10-15 % i.e. to 155-165 MeV, compatible with (i.e. below) QCD lattice values 
for the phase transition temperature from QGP to hadronic matter. Our value of 75 
is above 1 and this might look surprising. However, its value is decreased by some 15 
% when looking at a more realistic equation of state. In addition, we have imposed 
strangeness neutrality, it is possible that this is a too strong constraint when analyzing 
data taken in a very restricted rapidity region (see jEH] where a similar problem was 
encountered). Note that using a larger value of 7^, the size of the allowed window for 
initial conditions in flgure 1 increases. Aside of the uncertainty in the equation of state, 
there are other factors that influence the precise location and size of the window but 
few: value of the cross section, = 2 /m^, (taken constant for simplicity here) in the 
Glauber formula, value of the cutoff V = 0.5 which in fact is equivalent to a change in 
the cross section [2^1, and value of the initial time for the hydrodynamical evolution for 
which the canonical value of tq = 1 fm was assumed. 

With the initial conditions determined by flgure 1, particle spectra can be computed 
and compared with (all rapidity) NA35 data. This is shown in flgure 2. The agreement 
is reasonable. No decays have been included, this should in particular lead to an 
improvement of the low m± pion spectrum. 

Therefore there exist initial conditions of the hydrodynamical expansion such that 
NA35 strange and non-strange particle abundances and spectra can be reproduced 
simultaneously without extra assumption, in contrast to some of the freeze out models 
mentioned above. In fact it is puzzling that the continuous model leads to so many 
more pions than the freeze out model and it is necessary to investigate the reason. 

To illustrate more precisely the difference between (simple) freeze out and 
continuous emission scenarios, we compare in table 1 results from both scenarios, with 
To = 185MeV and fibo = 100Me\^ for the continuous emission case, Tf ^ut = 185MeV^ 
and fib f.out = lOOMeV" for the freeze out case (with the scaling factor Tjout in (2) taken 
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equal to 1 fm) ^. For both cases, we assume that initially the matter is a hadron gas and 
7s = 1.3. Heavy particles, in the continuous emission case, due to thermal suppression, 
are mostly emitted early |2ZI, i-e. in similar conditions than in the freeze out model in 
the case /if,o = and for some longer time if /x^o 7^ 0. Pions in the freeze out case are 
too few as discussed previously: 15.7 instead of 26 ± 1. Pions in the continuous emission 
case, on the other side, are emitted early and then on, and we get substantially more of 
them, 27 in agreement with data^. 

In a hydrodynamical model without shocks and dissipation, entropy is conserved. 
In the usual freeze out scenario, this is an important point because the initial entropy can 
be determined from the final multiplicity. To illustrate this connection, let us consider a 
thermalized massless pion fluid. In this case, the entropy density is related to the pion 
density at all temperature by s = 3.6n.,^. Therefore knowing the pion number at freeze 
out (from data) A^/'°"*, one can infer the entropy at freeze out S"-^'""* and the initial value 
of the entropy 

So = = 3.6Ar/-°"* (6) 

(since entropy is conserved). 

In the continuous emission case, we can compute the number of free particles in 
principle from (4) 

iVf- = Jd'xJ d'pD,{P^ffree) (7) 



d'^xD^ 



9 r,3P' 1 'P 



d'p'- 



{2TxfJ ^ £; e("-P)/^ - 1 1 - P 
rf^xa(x)n™=°(a;) 

where we used the factorization property arising from the zero mass assump- 
tion |2ni) and introduced a kind of weight, a{x) = D^{A^n^^)/n^^ with 
A^' = (/ dn/{47r) V/{1-V)J dVt sin 9 cos 0/(47r) V/{l-V)J dQ sin 6 sin 0/(47r) V/{1 - 
V)JdVt cos^/(47r)P/(l-P)). 

We re- write (7) by assuming that the fluid is divided in small volumes Va(r) moving 
with it and time is discretized with Tj+i = + Ar, so 

oo oo 

Nt'^ = 1111 ^r^V^{n)a{n)n2=\n) (8) 

a=l i=l 

At time Tj, AV"ci(rj) contains a mixture of free plus thermalized particles and we 
suppose that the thermalized pions break into those which have just done their last 
collisions (joining the free component of the fluid) and those which remain interacting. 
This interacting component is supposed to have reached thermal equilibrium at Tj+i 
and the process of separation between free and interacting part repeats itself. We note 
that both the separation (free/interacting) process and the re-thermalization one (even 

^ For both models, the values chosen for the parameters are typical. 

+ To compute the number of pions in the continuous emission case, the effect of continuous emission 
on the hydrodynamical evolution of the fluid is taken into account as in ,22| . 
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if incomplete) lead to entropy increase, which can manifests itself as pions. Exactly 
how much entropy is created and how much appears as pions is model-dependent, but 
this effect should always be present. Depending on the process involved, this can be 
a substantial effect. Indeed as shown in table 1 and discussed above, in the particular 
case of continuous emission with fint ~ fth (complete rethermalization), the amount of 
extra pions can even reach 70%. It would be worth checking if this is the case with the 
improved description of continuous emission [30]. The hypothesis fint ~ fth corresponds 
to the upper limit of entropy production. It is interesting to note that the observed pion 
abundance is reproduced just by the maximal entropy production within the present 
mechanism. Comparing results from UrQMD and the thermal statistical model jS], it 
was shown that the entropy density per baryon density increases by 15 % with increasing 
times, or alternatively for temperatures decreasing from 161 to 126 MeV, in a central 
cell 5 X 5 X 5fm^ for Pb+Pb at SPS. In our case, for such a central cell, away from the 
border, we have initially little continuous emission and the entropy per baryon increases 
slowly. The large increase of pion abundance obtained above occurs in the space-time 
domain where continuous emission becomes considerable. 

5. Conclusion 

In this paper, we discussed data on strange and non-strange particles at ultrarelativistic 
energy with light projectiles, from a hydrodynamical point of view. Some versions of 
the standard model with sudden freeze out can reproduce the strange particle data but 
underpredicts the pion abundance, if no extra assumption is made. In addition, it is 
usually necessary to assume two different freezeouts, chemical and thermal, to account 
for both strange particle abundances and particle transverse mass spectra. We showed 
that a hydrodynamical model with a more precise emission process, continuous emission, 
can reproduce both the strange and non-strange particle abundances without extra 
assumption in addition to being consistent with other types of experimental data such 
as transverse mass spectra. In the two freeze outs case, typical parameters are chemical 
freeze out temperature, baryonic potential and strangeness saturation factor as well as 
thermal freeze out temperature and baryonic potential; they are fixed by data. In the 
continuous emission, the parameters are: the initial conditions To, fibo and 7^, which are 
fixed by data; the average interacting cross section and initial time were chosen to have 
the canonical values 2 /m^ and 1 fm. We note that while in freeze out models, data 
give information only on the freeze out conditions, in continuous emission, observables 
depend on the whole fluid history. Therefore what observables teach us depend on the 
emission model. This point is reinforced by a comparison of Bose-Einstein correlations 
for freeze out and continuous emission pij 

Our main point is the following: in the usual freeze out scenario, a large pion 
number may be associated with a large entropy (cf. (6)). Here we showed that a 
large pion number can be generated by continuous emission. The reason for this 
increase is entropy generation during the separation and rethermalization processes 
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occurring during continuous emission. We stressed that how large is this increase is 
model-dependent, but in any case this possibility sheds a new light on the problem of 
pion emission at SPS. For freeze out, a large experimental value of N.^ implies a large 
initial entropy Sq and may be considered a hint of QGP formation (see e.g. [IHIISl)- 
For continuous emission, a large pion number and entropy may be a natural outcome 
independently of QGP formation. Consequently, a better understanding of particle 
emission in the hydrodynamical regime is necessary to assess the possibility of QGP 
formation in relativistic heavy ion collisions. 
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Table 1. Comparison of experimental particle abundances with continuous emission 
and freeze out predictions for S+S collisions at midrapidity. 
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Figure 2. Transverse mass spectrum computed for the initial conditions obtained in 
figure 1 and comparison with NA35 (all-y) data. No decays included. 



